Non-linear Response of the trap model in the aging regime 
Exact results in the strong disorder limit 
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I. INTRODUCTION 
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We study the dynamics of the one dimensional disordered trap model presenting a broad distribu- 
tion of trapping times p(r) ~ l/r 1+M , when an external force is applied from the very beginning at 
t = 0, or only after a waiting time t w , in the linear as well as in the non-linear response regime. Using 
a real-space renormalization procedure that becomes exact in the limit of strong disorder n — > 0, 
we obtain explicit results for many observables, such as the diffusion front, the mean position, the 
thermal width, the localization parameters and the two-particle correlation function. We discuss 

£C) ' in details the various regimes that exist for the average position in terms of the two times and the 

f*"^ ' external field. 
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Trap models provide a simple phenomenological mechanism for aging 0, 0, 0] . Their aging properties have thus 
\Q • been much studied, either in the mean field version 0, OJ El 0] > where "usual aging" occurs, or in the one-dimensional 
version 0,0, where both aging and subaging behaviors appear in different correlation functions. The mathematicians 
have also been interested by these trap models [ljj with special attention for the cases d = 1 0, ^| and d = 2 
[l3| . The one-dimensional version is moreover interesting on its own, since it appears in various physical applications 
concerning for instance transport properties in disordered chains [14l Il5| or the dynamics of denaturation bubbles in 
random DNA sequences |l6[ . 

The study of the response to an external field and its relation with the thermal fluctuations has been for some 
\& • years a central question in the description of the aging dynamics of glassy systems 0, ■ ^ i s thus natural to 
consider the trap models from this point of view. The studies on the violation of the fluctuation-dissipation relation 
in mean-field trap models have shown that the results depend on the observable 0] , and on the choice of functional 
form of the hopping rates |2fl l2l|. For trap models on an hypercubic lattice, there are no such ambiguities in the 
choice of observables and external fields, since the natural observable is the position : one is interested into the 
response of the position to an external bias, and in the thermal fluctuations of the position. Recen tly, the response of 
the one-dimensional trap model was studied via scaling arguments and numerical simulations in [22j, where various 
■ regimes were found depending on the relative values of the two times considered (t w ,t w +t) and the external applied 
field /, the main results being that in the linear response regime, the Fluctuation-Dissipation Relation (or Einstein's 
relation) is still valid in the aging sector, whereas the response always become non-linear at long times. We have 
' shown in |23| that these two response properties could be understood as consequences of a "non-linear Fluctuation 
^ , Theorem" that has for origin a very special dynamical property of the trap model. 

In this paper, we consider again the response properties of the trap model, but with a complementary point of view 
: we use a real-space renormalization group (RSRG) procedure to derive various explicit exact results in the limit 
of high disorder. The RSRG methods, that have appeared in the field of disordered quantum spin chains [24L |25|. 
have then been very powerfull to study the Sinai model [2(j, as well as reaction-diffusion processes in a Brownian 
potential |27| . In particular, since the out-of-equilibrium dynamics of the random field Ising model (RFIM) can be 
describded as a reaction-diffusion process in a Brownian potential for the domain walls, the RSRG method has been 
used 28] to study the response of the RFIM to an applied external magnetic field. The RSRG approach is also very 
appropriate to study one-dimensional trap models with a broad distribution of trapping times p{r) ~ l/r 1+M in the 
limit of high disorder /i — * 0, as explained in details in p9j for the directed trap model, and in [!j] for the symmetric 
trap model (i.e. in the absence of an external bias ) : the RSRG method is able to reproduce the exact exponents 
of the whole aging phase < fx < 1 and moreover allows to compute exact scaling functions for all observables in 
a systematic perturbation expansion in /x 0,H^]. In contrast with other usual methods for disordered systems, the 
disorder average is not performed at the beginning but at the very end : the RSRG procedure is defined sample by 
sample, all observables are then evaluated in terms of the relevant properties of a given sample, and can be then 
averaged with the appropriate measure over the samples. The RSRG approach thus provides a very clear insight into 
the important dynamical processes. 

In this paper, we generalize the RSRG approach describded in 9] for the unbiased trap model to include the 
influence of an external bias, and we obtain exact results for various observables in the high disorder limit /i — > 0. 

The paper is organized as follows. In Sectional we describe the trap model in an external force field /. In Section 
IIIII we explain the real space renormalization procedure in the presence of an external field. In Section llVl we describe 
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the effective dynamics when the external field is applied from the initial time t = 0, and we give explicit results for 
one-time observables in Section [V] In Section IVT1 we discuss the effective dynamics when the external field is applied 
only after a waiting time t w , and we compute the corresponding two-time observables in Section I VIII In Section 
IVIIII we discuss the various regimes existing for the disordered averaged mean position in terms of the two times and 
the external field. In Section iPXl we discuss the rare events that are responsible for the response, in the time sector 
where the effective dynamics gives no contribution. Finally, in Section El we compare the response and the thermal 
fluctuations in a given sample, to discuss the validity of the Fluctuation-Dissipation Relation. The conclusions are 
given in Section IXll and the Appendices contain more technical details. 



II. MODELS AND NOTATIONS 
A. Master Equation in a sample 

To study the trap model in an external force field /, we conside the following Master Equation 0, 



dt 



P">{x + l)W^ x} + Prix - - P$»(x) W«l x+1} + W[il x _ 1} (1) 



( f) 

with the initial condition P£i.q(x) = S Xt o- The hopping rates 



satisfy the detailed balance condition 



where the total energy 



w gl x±i} =e -^±H (2) 



U x = -E x - fx (4) 



contains both the random energy (—E x ) of the trap x and the potential energy {—fx) linear in the position x induced 
by the external applied field /. 

B. Law for the disorder 

The trap energies {E x } are quenched random variables distributed exponentially 

P (E) = 6{E > 0)^e-^ (5) 

This corresponds for the mean trapping time r = e@ E to the algebraic law 

q (r) = e(r>l)-^ (6) 

with the temperature-dependent exponent [i = ^- At low temperatures fi < 1, the mean trapping time J drrg(r) is 
infinite and this directly leads to aging effects. 

C. Link with the 'trap model with asymetry' 

To make the link with the 'trap model with asymetry h' studied in , we note that in the new time t = (2 cosh 
the Master equation (JIJ becomes 



d ItlH = q -V±P«\x + 1) + q -±W P P(x - 1) - 



dt t„+i * r, ■ ' 



(7) 



n-l 
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where 

9±{f)= e+ Pi +e -Pi = —2— (8) 

are the probabilities to jump on the right and on the left when escaping a trap, with the normalization q+ + q_ = 1. 
The asymetry 

M/) = <Z+-5-=tanh/3^ (9) 

varries between h(f = 0) = for the unbiased case and h(f — ► oo) — * 1 for the fully directed case. 

As soon as ft. > 0, the random walk is expected to become asymtotically directed on large scales. In this article, we 
will be interested into the case where the crossover towards the directed regime happens on large length scales, i.e. 
the local asymetry is very small h <C 1 or equivalently j3f <C 1. In this regime also considered in |23, the relation 
between the force and the asymetry is at lowest order simply linear 

h(f) = ^+0(W) 3 ) (10) 
and thus the results of the present paper can be straighfordwardly compared with j2^| . 



D. Entropy and Generalized free-energy 

As in similar models [23, El HH i the Shannon entropy 

S(t) = -J2Pt(x)lnP t (x) (11) 

X 

and the energy U(t) = J2 X Pt(x)U x (@J allow to define a generalized free-energy 

F(t) = U(t) - TS{t) =5^P t (s) [U x + TlnP t (x)\ (12) 

X 

The detailed balance condition implies that it is a non-increasing function dF }p < The equality with zero is 
possible only if all currents exactly vanish, corresponding to equilibrium. Here, since we consider the infinite line, the 
equilibrium cannot be reached and the free-energy will decrease with no bounds. 



III. REAL SPACE RENORMALIZATION PROCEDURE IN THE PRESENCE OF A FIELD 



We have already presented the real space renormalization procedure for the unbiased trap model in |9j. Here we 
wish to generalize this approach to the presence of an external bias / > 0. 



A. Notion of renormalized landscape at a scale R 

The basic idea of the Real Space Renormalization procedure [3, [53 |2^| is that the dynamics at large time is 
dominated by the statistical properties of the large trapping times. The renormalized landscape at scale R is defined 
as follows : all traps with trapping time r„ < R are decimated and replaced by a "flat landscape" , whereas all traps 
with waiting time r n > R remain unchanged. At large scale i?, the distribution of the distance I between two traps 
of the renormalized landscape at scale R takes the scaling form 

Mi) - 4^ (a = -L) (13) 



Rv V Rt " , 

where the scaling distribution is simply exponential V(X) = e~ x 

The distribution of the trapping times of the traps in the renormalized landscape at scale R is simply 

qn(T)=6(T>R)£(~Y (14) 

To relate the renormalization scale R to the time t, we have to study the time needed to escape from a renormalized 
trap. 
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FIG. 1: Definition of the escape time from a trap in the renormalized landscape : the trap of escape time to existing in the 
renormalized landscape at scale R is surrounded by two renormalized traps that are at distances 1+ and 2_ on each side. The 
escape time To is the mean time needed to reach either r+ or r_ when starting at tq. 



B. 'Escape time' from a renormalized trap to another renormalized trap 

We now study the 'escape time' T from a trap To existing in the renormalized landscape at scale R in the presence 
of a field / > 0. This trap is surrounded by two renormalized traps that are at distances 1+ and Z_ on each side (see 
Figure^). Whenever the particle escapes from the trap To, it escapes on the right with probability q+ = (1 + h(f))/2 
and on the right with probability g_ = (1 — h(f))/2 ©. 

If it escapes on the right, it will succeed to reach the trap r + with probability (|A10|1 

1 - e'Pf 

Vtif,l + )= 1 _ e . m+ (15) 

which varries between p^(l, f — ► 0) = 1/1+ for the unbiased case and p^(l, / — > oo) = 1 for the directed case. 
If it escapes on the left, it will succeed to reach the trap t_ with probability l|A13j) 

CI _ e -0f) p -(}f(i--i) 

pALI-)= [ - — Y^mi (16) 

which varries between p~ (/_,/—> 0) = for the unbiased case and pj(Z_, / — > oo) = for the directed case. 
Otherwise, it will be re-absorbed again by the trap To. So the total probability to escape when exiting from To reads 

p e {f, Z+, Z_) = —j —j (17) 



1. Number of sojourns in a renormalized trap before escape to a neighbor renormalized trap 

As a consequence, the probability En(n) to need (1 + n) successive sojourns in the trap To before the particle 
succeeds to escape either to the trap t_ or to the trap t + reads 

E R {n) =[l-p e {f,l + ,l.)} n p e (f,l + ,l^) (18) 

For large R, since we have l± = R^\± (fH5|l we obtain using l|AlllA14j) the scaling form for the probability of escape 

Pe(.f, i+,i-) * = PfR"; A+ = j±r, A" = l —) (19) 
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with the scaling function 



^ + ^>= 2(1 -V»+ )(1 - e -A-) <*» 



The expansion near F — * gives the first correction with respect to the symmetric case studied in 9] 



0(F;A+,A-) = 1 



1 f 

A+ + ~ 



X++X ^-^ 3 r^ 3 F^O(F-) (21) 



24 1440 



whereas in the other limit F — > oo, we have 

F 



0(F; A+, A") ^ £■ (l + e-^ + + e"^ + ...) (22) 



As a consequence, at large scale R, the number n of returns is distributed exponentially 

1 

<n> (R, F) 

where the mean number of returns 



E R (n) ~ 1 e ~<">w^) (23) 



varries between 



and 



<n> < fl ' F ' = «i?^F) (24) 

o\+\ — 

<n >(R,F^Q) = W^ rTJ - (25) 

<n>(i2,F^ 0O )~^| ; = A (26) 

2. Total time spent inside a renormalized trap before escape to a neighbor renormalized trap 

Let us now consider the probability distribution Pi n {ti n ) of the total time tj n spent inside the trap To before its 
escape. It can be decomposed into the number n of sojourns, where n is distributed with l(T%|l 

l+n 



tin — ^ t n (^^) 



where t n is the time spent during the sojourn i in the trap To, so it is distributed with the exponential distribution 
with mean time To- Actually, since n is large in the large R limit, we have the central-limit theorem 

tin — n <tn >— nTQ (28) 

n — >oo 

Since the number n is distributed with exponentially (|23|1 . we finally obtain that U n is also exponentially distributed 

tin 

P in (t in ) ~ Tv (29) 

with the characteristic time 

T = t < n > (R, F) (30) 

Since the smallest trapping times existing in the renormalized landscape at scale R is To = R, the time spent inside 
the trap To before it succeeds to escape scales as 

tMF) RZoo R<n>{R,F) = (31) 
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3. Total time spent during the unsuccessful excursions before the escape 

Among the n unsuccessful excursions, there are to excursions on the left and (n — to) excursions on the left, where m 
is distributed with the binomial distribution 2 _Tl C™. Since n and to are large, we again have a central-limit theorem 

m n—m 

tout = £ *» + E ^ ~ m0 us(f, 1-) + (" " m)9+ s (f, 1+) (32) 

n=l j=l 

where 0^ s (f, I) represents the mean time needed to return to when starting at (±1) without touching the point (±Z) 
in a fiat landscape. Using the asymptotic behavior (|A17IA21|) 

9+ s (f,l) ~ 10(u = /?/!) (33) 

I— >OC 

with the scaling function given in (JA18J) . the ratio between t ou t and U n (|28[1 scales as 

^(i?, F) ~ ^ s(/ ' L) + ^ s(/,?+) ~ fl"" 1 [6(FA+) + 9(FA-)1 (34) 

tin 2tq 



which varries between (|A19(I 



and 



— (R,F^0)~ RV- 1 (35) 



t -^i { R,F-,oo)^R^^^- (36) 



We may thus neglect t out with respect to t m at large scale R. 



4-. Time spent during the successful excursion to escape 

We finally consider the diffusion time i^r* of the successful escape to the neighbor renormalized landscape (±Z) 
when starting at x = 1 without visiting x = 0. The mean time needed to reach x = l + when starting at x = 1 for a 
random walk conditioned not to visit a; = takes the scaling form l)A22l 

t+ ff (f,l + )~l 2 + D(u = 0fl+)) (37) 
where the scaling function D{u) is given in l|A23|) . Using also (|A26|) . we obtain the scaling 

t -^(R, F) ~ R^ViFX) (38) 

tin 

which is the same as (15*4^1 and thus t<n //is again negligible with respect to t in (|3*TJl . 



5. Conclusion 



So we obtain that the total time 

tesc tin ~t~ tout tdiff (39) 

needed to escape is actually simply given by the time ti n spent inside the trap to . So the distribution of t eS c is given 
by the exponential l|29|l with the escape time Jo- 
in conclusion, a trap of the renormalized landscape at scale R has a trapping time r distributed with (|f4|) , but has 
an 'escape time' proportional to r 

T = r<n>(R,F)=r w ^ ) (40) 
in terms of the two rescaled distances X± to the neighbors. 
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C. Distribution of 'escape times' in the renormalized landscape 



The distribution of the escape time T in the renormalized landscape at scale R is obtained by averaging over r, A± 

+00 / n\ f r+OO rp/J 

*(T;f) - I J -»-,(!■- j^^) (41) 



It thus takes the scaling form 

Qr(T; f) =^Q,(f=^;F) (42) 

where the scaling function 
presents the asymptotic behavior 

Qu(T-F) ~ , M c M (F) (44) 



(43) 



with the prefactor 



/■ + OO _ _ /- + 00 

c„(F)= d\+d\- e - x - x [cj)(F;\ + ,\-)] fl = l- f i d\+d\- e - x - x \n[<p{F;\+,\-)] +<D(n 2 ) (45) 

Jo Jo 



For small fi, the probability distribution ^T, _Fj is dominated by its long tail, and we may approximate it by 

Q, (f ) ~ 9(f > f,{F))^ (^jP-^J (46) 

The cut-off T M (F) chosen to preserve the normalization is determined by the coefficient in the long tail part (|44I45|I 

f^F) = (c^F))i =e s ^)+0( M ) (47) 
where the function S(F) obtained from l|45|) can be computed with the explicit expression 1|2U|) 

r+°° + 

S(F) =- d\+d\- e - x - x ln[<f>(F;\ + ,\-)] (48) 







lnF + ln2- 7E -V ( 1 + J - I 1 + i ) ( 49 ) 



F J F \ F 

in terms of ip(z) = T'(z)/T(z) and the Euler constant 7^ = 

For the unrescaled probability distribution l|42|). this corresponds to the cut-off 

T (i?,/) = i? 1+ ^e s(F) (50) 

D. Choice of the renormalization scale i? as a function of time 

The renormalization scale R has to be choosen as a function of time by the requirement that this effective cut-off 
Xo(i?, /) is exactly t, meaning that at time t, only traps with escape times T > t have been kept, whereas all traps 
with escape times T < t have been removed and replaced by a flat landscape. So given a renormalisation scale R and 
an external field /, the corresponding time t of the dynamics reads 

t = t(RJ) = R 1+ »e s ( F =M R ^ (51) 
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The renormalization scale R(t, f) as a function of time t and field / is thus defined by the implicit equation 

t = R 1+ »(tJ)e s W R " (52) 
At short times, assuming F — (3JR^(t) < 1, we may use the expansion 

1 



S{F) = In 2 - 1 - j E - j^F 2 + 0(F 4 ) 
to solve (|52() to obtain the first correction with respect to the unbiased case 



R(tJ) = 



t \ 1+ " 



12(1 + m) \T 0j 



where Tq — e s ^ — 2e 1 7E . This solution is thus valid for 

where the time scale t^{f) was already shown to play an important role via other approaches |22l |23| ]. 
At long times, assuming F = PfR^(t) ^> 1, we may use the asymptotic behavior 

W = -lnf-^+0(i,) 



(53) 



(54) 



(55) 



to solve <|52|) which leads to 



R(tJ) = ^-t 



(56) 



(57) 



which is valid for t ^> t fl (f) (|55|l . 

As a consequence, the characteristic length scale corresponding to the mean distance between renormalized traps 
at scale R(t) (fT^jb 



e(t,/) = fl"(t,/) 



behaves as in the unbiased case at short times 



and as 



at long times. 



£(*,/) * t— [l + 0(/x)] 



at j) * jpftrii+oQj)] 
*>t^(/) 



(58) 
(59) 

(60) 



IV. EFFECTIVE DYNAMICS IN THE LIMIT /x -> 
A. Probabilities to escape on the right or on the left in the renormalized landscape 



We are now interested into the relative probability to escape on the right rather than to the left which reads (|15I17(I 
using l|AHIA14jl 



u>+ (/,!+,*_) = 6+ ^ (/ ;V - W + (F, A+ , A- ) 

Pe{J, l+,l-) 



with the scaling function 



W+(F,\+,\-) 



1-e 



-FA" 



1 _ e -F(A+ + A-) 



(61) 



(62) 
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The complementary probability to escape on the left reads 

e - FX ~ (I - e - FX+ ) 

W-{F,X + ,X-)= 1 „ 6 Ua ++ a-) J (63) 
The expansion near F — ► gives the first correction with respect to the symmetric case Q 

F A+ - A" 



W+(F,\+,\-) = A 



A- + A+ 



1 1 

A+ + ~ 



2 12 



1 1 

A+ + ~ 



F z + 0{F A ) (64) 



whereas for large F, the first corrections with respect to the directed case are given by 

W+(FX+,\-) ~ l- e - FA "+e- F ( A++A ^+... (65) 

F— >oo 

B. Rules for the effective dynamics 

We thus define the effective dynamics by the following rules : 

The particle starting at the origin O will be at time t either in the first trap M+ of the renormalized landscape at 
scale R(t, /) on its right or in the first trap M_ of the renormalized landscape on its left. The weights of the traps 
M + and M_ are given in terms of 162I63|I by 

P[m_m + ](M+\0) =W + (F,X+,\-) (66) 

P[m_m +] (M_|0) =W_(F,A+A-) (67) 

We now verify that this effective dynamics presents some important properties. 

C. Consistency upon iteration 

The rule for the effective dynamics is consistent upon iteration. Indeed, suppose there are three consecutive traps 
: the trap M_ is at a distance i_ from the origin on the left, the trap M + is at a distance 1+ from the origin on the 
right, and the trap M ++ is at a distance I from the trap M + on the right. 

Suppose that the trap M+ is decimated before the traps M_ and M++. The new weights for the traps M_ and 
M ++ become 

Pm. =P[m_m + ](M-\0)+p [m _m ++ ](M-\M + ) P[m _ m+] (M + \0) =p [M _ M++ ](M-\0) (68) 

1 _ e -Fi- 

Pm++ = P[M-M ++ }(M ++ \M + )p [M _ M+] (M + \0) = - F{l++l _ +l) = P[M-M ++ ]{M ++ \0) (69) 

and thus the rules l|67(l for the occupancies of renormalized traps are consistent upon decimation of traps in the 
renormalized landscape. 

D. Conservation of the mean position in any sample for / = 



As already emphasized in |22j, |23j , the trap model has a very special property : when there is no external field / = 
the mean position is a conserved quantity in any given sample, as a simple consequence of the master equation Q 
for any realization of the trapping times {t„}. Here in the effective dynamics, the mean position vanishes indeed in 
each sample with the weights in zero- field / = given by (|64|l . 



E. Non-linear fluctuation theorem 

The non-linear fluctuation theorem obtained in |23| says that, in any given sample of the trap model the diffusion 
front Pj (n) in the presence of the external field (+/) and the diffusion front (n) in the presence of (— /) are 

related by the extremely simple property 

y \ } n) = (70) 
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In the effective dynamics, this property is satisfied since the renormalization landscape at time t is the same for (+/) 
and (— /), so the two possible positions M + and M_ are the same, and the weights (|62I63[) satisfy the properties 

W+(-F,\+,\-) =e- FX +W+(F,\+,\-) (71) 
W-(-F,\ + ,\-) =e FX -W-(F,\+,\-) (72) 



V. OBSERVABLES WHEN THE FIELD IS APPLIED FROM THE BEGINNING t > 

A. Diffusion front in a given sample 

In the effective dynamics, the diffusion front in a given sample takes the scaling form 

P^in-J) = -L^°> (X = 1 £ r) ;F = pmtj)) (73) 

where the characteristic length scale /) has been defined in (|52I58() from the renormalization landscape, and where 
the scaling function reads according to the effective dynamics (|67|) 

-F\-(i_ p -F\+\ 1 _ -FX' 

in terms of the two rescaled distances X± between the origin and the nearest renormalized traps. The joint distribution 
of the two rescaled distances is completely factorized 

X>(A+,A-) = 9(\+)0(\-)e- x+ - x ~ (75) 

We now compute the various observables that can be obtained from the sample-dependent diffusion front (|74|l and 
the measure l|75|) over the samples. 



B. Disorder averaged diffusion front 



The averaged diffusion front takes the scaling form 

PM) = j^g, (X = ^ ,F = m t, /)) (76) 
where the scaling function reads at lowest order in fi — * 

go(X,F) = e— 

= e -\x\ 

which interpolates between the diffusion fronts for the unbiased case |9( as F — > and for the directed case |29| as 
F — > oo. Here, as for other results below, we have given two equivalent forms for the diffusion front : the first one 
with the integral is more appropriate to study the small F behavior, whereas the second one with the infinite series 
is more suited for large F. 
The simple property 

go(-X,F) = e- FX g (X 7 F) (79) 



0(X > 0) + 6(X < 0)e- F|x| | / d\e~ x - 

- 1 Jo 1 

9(X > 0) + 6(X < 0)e- F l*l] Jl 



1 - e 



FA 



e -F(|X|+A) 



> (l + Fn)(l+F + Fn) 



-nF\X\ 



(77) 
(78) 



is actually expected to be valid for arbitrary fi as a consequence of a non-linear Fluctuation Theorem, as discussed in 
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C. Disorder averaged mean position 

The disorder average of the mean position takes the scaling form 

+00 

<x(tj)>= J2 xP^J)=t{t,f)X^F = pfZ{t,f)) 

X——00 

where the scaling function X^ reads at lowest order in fi — > 



r+oc r+oc fx cotll EA _ 1 

X (F) = I dXXg (X,F) = / dXXe " 



1 1 



-A 2 — 2 1 = l_±_ _Lw/' + - 

F F F 3 \ F 



The series expansion for small F reads 



F F 3 F 5 , 7n 



whereas the asymptotic behavior for large F reads 

+00 



^ (1 + Fn) 3 

n— 1 v * 

Using _F = /3/£(t, f) and l|59l6Uf) . we finally obtain the expansions 



F F 3 15F 4 \F 5 



0f ( t v+* 



<*(*,/)> = ^ 
at short times, i.e. for i <C ^(/) (|5 5|) . and 



2m 

1~\W? (^) 1+ %0((/3/t^) 4 ) 



<*(*,/) > 



o 



pf v(/3/) 2+ ^) 



for t ^> tfj,(f). The leading terms are thus in agreement with the scaling analysis presented in [2^ 



(80) 



(81) 



(82) 



(83) 



(84) 



(85) 



D. Thermal width 

The disorder averaged thermal width l|74|) presents the scaling form 



+00 



<Ax*(t,f)>^ x2p t&f)- 

X— — OO 

where the scaling function reads at lowest order in [i 



+00 



-1 2 



^e(tj)A^(F=pmtjy) 



(86) 



A (F) = 



+ OO 



d\e- x 



AH 



coth ■ 



FA 



FX 
2 



2 sinh 2 Q 



F F 4 V F ) F° \ F 



(87) 



It is thus directly related to the scaling function Xq(F) governing the mean displacement ijSJ via the simple relation 



A (F) = 



1 A 

F + dF 



X (F) 



Using the asymptotic behaviors for small and large F 



A (F) = 1 - -F 2 + F 4 + 0(F 6 ) 
Am _ 1 4C(3) 7T 4 / 1 



(89) 
(90) 
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with F = f3f£(t, /) (|59l 160(1 . we obtain the following leading terms for the unrescaled thermal width 



<Axi(t,f)> = 



. 2 m r 
t \ 1+ » 



(91) 



for t < t^f) ijSSJ, and 



<Ax*(t,f) > 



P 1 



4C(3) 



O 



(92) 



for t > t M (/). 



E. Mean square displacement 



In addition to the thermal width (JS5J, it is also interesting to consider separetely the two terms, namely the mean 
square displacement 



2 2 



Po(F)^<^> = 2-^ + ^ + -^[l + -l=l , ? 



1 



F 2 F 4 



0(F 6 ) 



and the sqare of the mean displacement 



F F 2 F° \ F ) 3 



0(F 6 



(93) 



(94) 



in the rescaled variable X = x/£(t, /) as usual. 



F. Localization parameters 

The localization parameters, defined as the disorder averages of the probabilities to find k independent particle in 
the same trap, follow the scaling form 



X — — QO 



where the scaling function at lowest order in fi reads (JT^J) 

"+OC 



yi 0) (F) 



ii Jo 



-FX~ 



I _ e -F(X++X- 



1 - e 



-FX~ 



I _ e -F(X++X~) 



+ 00 



^ fcr (± + m)T(k + m) 
^ FT(1 + m)r (1 + i + k + m]) 



1 + Fm l + F(k + m)_ 
The series expansion for small F gives the first correction with respect to the unbiased case Q 



2k(k- 1) 



F 



2k(k- l)(fc 2 - 15k -A) 



k + 1 (k + l)(k + 2)(k + 'S) (k + l)(fc + 2)(fc + 3)(fc + 4)(fc + 5) 

whereas for large F, the first corrections with the directed case [2^ read 



4V) =l-(lEuler+m)j + 



1 



{lEuler+^{k)) 2 , 7T 2 3 



1 



4 2^J^ + °VF3 



(95) 



(96) 



+ 0(F 6 ) (97) 



(98) 
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For k = 2 and k = 3, we may moreover sum the series (|96|l to obtain 

y?\F) =i 
y<°'(F) i 



<">,n-, , 3 3 --L/fl + 1 



2F F 2 F 3 



F 



which are thus related by the very simple relation 



i-^ 0) (f) i-y^(F) 



(99) 
(100) 

(101) 



This relation actually reflects the the two-delta structure of the diffusion front in any given sample : indeed, if we 
note pi and p 2 the weights of the two delta peaks with pi + p 2 = 1, it is immediate to obtain the relation (|101fl since 



1 - Y 2 (0) = ^ + p 2 f - p\ - p 2 2 = 2 PlP2 

1 - ^ 3 (0) = (Pi +P2? -P\-Pl= 3 P1P2 (P1 +P2)= %PlP2 



(102) 



G. Entropy 

The entropy ljll|l is closely related to the localization parameters ()95(l since it reads 

= 5^(F = /?/£(*,/)) 



S(t,f) = - ]T Pt(x) In P t (x) 



dk 



where the scaling function reads at lowest order in \x 



5(0) (j?) 



F 2 



+ °° / 1 
? [[ (I + F-\ 



(l+F + Fm) 2 (1 + Fm)' 



if> (1 + to) — ?/> ( 1 + — + m 



The asymptotic behaviors are given by the series 



for small F and by 



i p2 p4 

5(°)(F) = + — + 0(F 6 ) 

K ' 2 12 20 K ' 



S (0) (F) -** ^Knf 1 



(103) 



(104) 



(105) 



(106) 



for large F. So for the unbiased case / = 0, the entropy remains forever frozen at the value S^(F = 0) = 1/2. For 
the biased case, the entropy decays towards zero as the directed character gets stronger as time grows. 



H. Two-particle correlation function 

The two-particle correlation function reads 



^oo +oo 



C(l,t)= J2 E W^'K|*-x'L- Y2 0) k0 + Jf^rC 



X — — oo x' — — oo 



at j) 



(107) 



where the weight of the delta peak corresponds as it should to the localization parameter discussed above, 
whereas the scaling function of the long-ranged part reads at lowest order 



Co (A, F) = 



— A 



1 - 2F\e- FX - e- 2FX 



F(l-e- FX ) 2 

In particular, the first correction with respect to the unbiased case reads @ 



4 0) (A,F) = e 



-A • 



2\2 



1 - 



F 2 X 
~3CT 



4 \4\ 



0(F 4 A 



(108) 



(109) 
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VI. EFFECTIVE DYNAMICS WHEN THE FIELD / IS APPLIED FOR t > t w 

We now study the dynamics in the following aging "experiment" already considered in |22j| : the system first evolves 
with no external field / — during the interval t G [0, t w ], and then an external field / > is applied for t S [t w , +oo[. 



A. Time sector governed by the effective dynamics 

The scale of the renormalized landscape corresponding to time t w and to no external field / = reads (|54|) 

i 

R(t w ,f = Q)= (HO) 
The corresponding length scale between renormalized traps is given by (|59|1 

e(w = o) = an) 

The state reached at t w , which is made out of two delta peaks (|74|l , has to be considered as an initial condition for 
the dynamics in the presence of / > in the new time (t — t w ). Since at t w , the particle is typically in a trap of trapping 
time t > R(t w ,f = 0), the effective dynamics corresponds to no move as long as R(t — t w ,f) < R(t w ,0), whereas 
the effective dynamics governed by the decimation procedure becomes "active" for R(t — t w ,f) > R(t w , f = 0). As a 
consequence, it is useful to introduce the parameter 

,,, ,v_ Z(t-t w ,f) R»(t-t w ,f) 

a(t ^ f) = «w = o) = mt w ,f = o) (112) 

that measures the ratio of the length scales of renormalized landscape at the two scales R(t — t w , /) and R(t w , / = 0). 
The response in the time sector a < 1, which is governed by rare events will be discussed separately in the section 
IIXI In this Section, we will concentrate on the time sector a(t, t w , /) > 1, where the effective dynamics governed by 
the decimation procedure is the leading effect. 



B. Role of the characteristic time t M (/) associted to the external field 

We have already seen that a bias / introduces the characteristic time scale = [jj) " f55|l - It is thus natural 

that the domain in (t — t w ) corresponding to the time sector a > 1 actually depends on the relative value of t w with 
respect to t^(f). 

For t w -C t^(f) 7 the time sector a > 1 corresponds to the domain (t — t w ) > t w , and the ratio a interpolates between 

a(t,t w ,f) =(^^y + " for t w <t-t w <t M (f) (113) 

and 

fii(t-t )V 

a(t,t w J) = 2 \ for t-t w »t M (f) (114) 

1 

On the other hand, for t w t^(f), the the time sector a > 1 corresponds to the domain (t — t w ) > -^jtw +>1 , and 
we have ()114f) in the whole time sector a > 1. 

To simplify the reading in the following, we will use the simplified notations 

R = R(t — t w , /) 
R w =R(t w ,f = 0) 

R^ L 

a = a(t,t w ,f) = — jr (115) 

lt W 
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C. Statistical properties of the renormalized landscape at two successive scales 

To study the two-time effective dynamics, we will need the probability measure for the renormalized landscape at 
the two successive scales R w and R. Since the two sides of the origin are independent, we first consider the half line 
x > alone. The joint probability that the first renormalized trap at scale R w is at distance Z+, has a trapping time 
t+ and that the first renormalized trap at scale R is at distance l + reads 



A Rw ,n{T+,lt,l + ) = G{<>R w )- 



i± 



(T+)!+^ 

After integration over the trapping time rj, we obtain the scaling form 

A R „, R (l+;l+) = J dT +A Rw , R (T+J+ ] l+) = 1 ^A 
with the scaling function 



r+ > R)S(l+ - 1+) + 9(t+ < R)6(l+ > l+)—e 



1 



7+ /+ _ /+ 

w ~ RH, ' - R» ' RZ, 



A(X+;X + :,a) = 



-6(X+) 



1 - - ] 6>(A + > 0)e" 



Since we have the same properties for the half- line x < 0, the measure for the full line reads 



(116) 



(117) 



(118) 



= 0(A+ > 0)9(\- > 0)e" 



1 - - ) 0(A+ > 0)e 



— A" 1 



-<J(A" 



(119) 



1 - - ) 6>(A~ > 0)e 



-A" 



D. Rules for the two-time effective dynamics 

Given the configuration (A - , A~, A+, A + ) charaterising the renormalized landscape at the two successive scales, the 
two-time diffusion front takes the scaling form 



with 



PyX) t , X w , t w 1 0, 0) 



1 .pW 



R&R* {>>-A»M 



(120) 



V 



(o) 

{A",A~,A+,A+} 

— 5(X W — A^) 



(X w ,Y,F,a) 



W+{\+, \- + X ™ +Xw ,F)5(Y - A+) + T4^_(A+, A" + A ™ +A ™ ; F)£(Y + A" + A ™ + A ™ ' 



A+ 



Am + A w 

VK+(A+ 



-s(x w + \-) 

w w , \-,F)5(Y - A+ - m ; 



W_ (A+ + Xw + Aw , A", F)5(y + A~) 
a 



(121) 



We are now in position to compute various observables in the aging regime from the knowledge of the sample- 
dependent two-time diffusion front l|122[l and from the measure (I12U|I over the samples. 
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VII. OBSERVABLES WHEN THE FIELD / IS APPLIED FOR t > t v 
A. Disorder averaged two-time diffusion front 

The disorder averaged diffusion front with the measure (|12U[) reads 



/*+oc /*+oo /*+oo P 

VW(X w ,Y,F,a)= d\+ d\- d\+ 

Jo Jo Jo Jo 



0(X W > 0)0(Y > 0) + e- F ^6(X w < Q)6(Y < 0) <?<°{ (\X W \, \Y\,F,a) 



d\-A (A+; A" ; A+; A" ; a) V™_ A _ ^ x+} (X w , Y, F, a) 

(122) 



0(X W < 0)9(Y > 0) + e- F ^9(X w > 0)0 (Y < 0) (\X W \, \Y\,F, a) 



,(°) 



where 



gf + {X w ,Y,F,a) =e 



-5{Y)+(l--\, ' 

a \ a 



+oo \ — 

i) X w + X. w 



1 \ 1 - e 
a) 

1 - e 



1 

1 

"/ Jo 



cZA^e^" 



1 - e 



1 - e V 



(123) 

(124) 
(125) 



and 



+oo \ + 



+ 00 

1 - — ) / d\- e - x ~ 



1 - e 



-FA" 



/o A TO + 

The form (|123|l presents the simple property 



a a 



1 - e -F(Y+\~) 



1--)9(Y 
a 1 



A w +X W _ (Y - 



a 



p(°> (--X*, -F, F, a) = e"* F, F, a) 



(126) 
(127) 

(128) 



which is expected to be true for arbitrary [i, as a consequence of the non-linear Fluctuation Theorem discussed in 



B. Law of the relative displacement Y between t w and t 



We now consider the partial law of the rescaled relative displacement Y between t w and t, i.e. we integrate the 
disorder averaged diffusion front 1)123(1 over the position X w at time t w : we obtain the form 



V(°)(Y,F,a)= / dX w V^ (X w , Y, F, a) = 6(Y > 0) + e~ F ^e(Y < 0) G (Y, F, a) 



where 



p + oo 

G (Y, F,a)= / dX u 
Jo 



gf + {\X W \, \Y\,F,a) + (\X w \,\Y\,F,a) 



(o) 



Using the intermediate results 



11 r+°° i _ p—Fu 

dX w gf + (X w , Y, F, a) = —5(Y) + - e - Y J du — — [ e "« - e""*] 



and 



+oo 



I dX w gM(\X w \,\Y\,F,a) = ^e- r -e~ ay )J due^^—^—^ 



i r+oo 
' -y „-aY\ 



1 - e 



-Fti 



(129) 



(130) 



(131) 



(132) 
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we finally obtain 



V^(Y,F,a) = -S(Y)+ 6(Y>0)+e- F ^9(Y < 0) G ns {\Y\,F,a) 



where the non-singular part is given in terms of the function 



G ns (Y, F, a) 



y r+°° 



du 



1 



1 - e -F{Y+u) 



(2 - e- {a - 1)Y )e- u 



+ 00 



-FnY 



H 



(1 + Fn)(l + F + Fn) 



(2 



H 



(a + Fn){a + F + Fn) 



(133) 

(134) 
(135) 



The presence of a singular part in 5(Y ) for the rescaled relative displacement was already found in the case of the 
Sinai model by the RSRG approach |26j and by mathematicians |35j. Here the weight of S(Y) is simply given by 
(1/a) l|112|) which represents the probability for a renormalized trap at scale £(/ = 0,t w ) to be still present in the 
rcnormalized landscape at the new scale £(/, i — t w ). 



C. Disorder averaged mean position 



Since < x w >= in any sample, in the effective dynamics in the limit \i — > but also more generally for any fi as 
a consequence of a special dynamical symmetry 0, the mean position < x > at time t in a sample also corresponds 
to the mean displacement < (x — x w ) > between the times t w and t. The disorder average of the mean position takes 
the scaling form the scaling form 



<x>(t, t w J) = RfX^F = 0fR», a = S) 



where the scaling function is simply given in terms of the law of the relative displacement Y (|133fl 



dv- 



1 - e 



-F'L 



F F 3 \ F 
1 

' (a - l)((a - l) 2 - F 2 ) 



dYYe- y (1 -e- rr )(l -e**-* v ) (2 - e 
111 1 



2a 2a 2 2(a — 1 + F) 2a 2 (a — 1 — F) 



(136) 

(137) 
(138) 
(139) 



We now consider various limit expressions. 

On one hand, for fixed a, the asymptotic expressions for small F and large F which generalize (|82I83|I read 



and 



X {F,a) 





F 


1 


1 " 




1 
















Hi 


a _ 


~2 


6 






" 6 





1 3 1 

5 V a a a° 



F 5 + 0(F 7 ) 



-A. ± -JLVifi-iU-i 



2a 2a 2 / 2 



a 2 F 



1 



1 



(140) 



(141) 



On the other hand, for fixed F, we have at the beginning of the aging regime where (a — 1) is small 



X (F,a) =(a-l) 



1 



1 



+ 0((a-l) 2 ) 



2 F 2F 4 

whereas for large a, the first corrections with respect to the previous result l|81|l corresponding to a — > oo read 

F 



*o(F,a)=l-l-^Wl+i 



2a 2 



F 2 ( 1 
-F + T+ ^' 1 + - 



4 + 



(142) 



(143) 
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D. Mean-sqare displacement during [t w ,t] 
Similarly, the mean-sqare displacement during [t w , t] takes the scaling form 



D{t, t w ) = < (x(t) - x{t w )) 2 > = R^V^F = 0fR», a = - fr ) 
where the scaling function is again obtained from the law of the relative displacement Y (|133fl 



V {F,a) 

2 2 



dYY 2 [l 

1 1 

F ' F 2 a ~ a 3 
4(a- 1) 



-FY1 



du 



1-e 



-Fv 



1 - e -F(Y+u) 



1 



(2-e 
1 



-(a-l)y-j e -u _ g-au 



-1 + F (a -1 + F) 2 a 2 (a-l-F) 2 a 3 (a-l-F) 



+ 



F((a - l) 2 - F 2 ) 
2(a — l) 2 - F 2 









-?)] 





(144) 

(145) 
(146) 
(147) 



F'Un -|)'i_F 2 )^ ( 1+ F) + F 2 ((a-1) 2 -F 2 ) 



*»( 1 + j) 



For fixed a, the asymptotic expressions for small F and large F respectively read 
1 



V (F,a) = l- 



112 1 

3 ~ 2q^ + 3q^ ~~ 2o^ 



F 2 + 



11 1 2 11 

3 + 6^ + 5^~5o :l+ 5Q : 5 + 6a6 



and 



A,(F,a)=2---i- 1~U 



2-a- 



1 \ 1 







1 

F 3 " 



a 3 / F 2 

For fixed F, we have at the beginning of the aging regime where (a — 1) 

(a-l) + 0((a-l) 2 ) 



F 4 + 0(F 6 ) (148) 



(149) 



2?o(F,a) 



4 3 2 1 ,////- 1 



whereas for large a, we obtain the first corrections to the previous result (|93|) 
£>o(F,a) = 



2 2 2 „ , 



-2 + F 2 + ^/ + JLtf' ( 1 + A 



F 



F 



F 2 



F 



4+°(A 

ar V a 



(150) 

(151) 
(152) 



E. Thermal width 



The rescaled thermal width 



A M (F,a) = < X 2 > - < X > 5 



(153) 



may be similarly computed, but since the full expression is rather lengthy, we will only give the asymptotic forms. 
For fixed a, the asymptotic expressions for small F and large F respectively read 



Ao(F,a) 



1 



2 1 1 

3 + ~ 3a 4 



F 



1 



and 



A (F,a) 



1 1 1 

3a 3a 2 3a 3 



2 2 \ 1 

3 ~ 3a 3 J F 



3a 4 9a 5 9a 6 



a 1 



F 4 + 0(F 6 



1 



6 6 6a 2 6a 3 / F 2 
For fixed F, at the beginning of the aging regime where (a — 1) is small, we have 



O 



1 

F 3 



(154) 



(155) 



Ao(F,a) =1 



2 2 2 , ... / 1 \ 1 , .... / 1 
-2 + -ib'" 1 + — t4> ' 1 + — 

F 3F 2 3F 5 V \ F 3F 6 ^ \ F 



(a-l) + C>((a-l) 2 ) (156) 
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whereas for large a, we obtain the first corrections to the previous result (|87() 

1\ 1 ..„/. 1\ F 2 



A (F,a) 



F 



F 5 \ F 



3 3^1 F 



For the case £ w = 0, corresponding to a = oo we have previously found a very simple relation l|88|) between the two 
scaling functions describing the mean position and the thermal width. Here, in the aging case t w > with finite a, 
there doesn't seem to exist a simple generalization of ((8811 . 

VIII. VARIOUS REGIMES FOR THE RESPONSE IN THE LIMIT /x -> 



In this Section, we translate our results for the disorder averaged mean position 1(1361139(1 into the original unrescaled 
quantities, using the definitions for F — flf£(t,f) and for a(t,t w ,f) 1)112(1 . We have to distinguish various regimes, 
according to the relative values of the times (t,t w ) and of the characteristic time scale t M (/) (|55jl associated to the 
force /. 



A. Case t w < t„(f) 



For t w <C t^(f), we obtain from I|14(JI113[I the behavior in the sector t w < t — t w <C ^(/) 



Bf m 

< x >{t,t w , f) ~ ^L(t - t w )T+t 

and from 1(1411114(1 the behavior at long times t — t w ^> t^{f) 



1 - 



t tn 



2m 
1+H 



< X >(t,t w , f) 



(t - t w y i 



2 M N 



(158) 



(159) 



Moreover, at the beginning of the effective dynamics regime a — > 1 (|113fl . the expressions (|140fl and H142|) coincide 
and give 



< x >(t,t w , /) t _ t -_^ 



whereas asymptotically when a — > oo HI 14(1 . the expressions l|141|l and l|143|) give 



(160) 



(161) 



B. Case i„ > i M (/) 



For in, ^> tf_t(f), the time sector a > 1 1(114(1 implies that t — t w ^> t fi (f) and we obtain from ((1411114(1 the behavior 

i 

at long times in time sector t — t w ^> and t — t w > -^jtw + ^ 



(162) 



Here at the beginning of the aging regime t — t w — > -Jrii, 4 ^ , the expressions ((141(1 and ((142(1 coincide and give 



< x >(t,t w , /) 



AT 

^ j- T+7T 

i 2 

— i— TT ^ 



- 1 



(163) 



whereas at at the end of the aging regime a — > oo ((114(1 . the expressions l(141|) and l(143|l again give ((161(1 . 
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C. Discussion 



The results of this Section, that rely on the effective dynamics picture, are valid at lowest order in the high disorder 
limit fi —* 0, and in the asymptotic aging regime where the two times are big t — > oo, t w — > oo, with the parameter 
ot(t,t w ,f) > 1 (I112|l being fixed. We have found various interesting behaviors depending on the relative values of 
the parameters (t,t w ,f). Whenever they can be compared, our results agree with the scaling analysis presented in 
|22|. In the next Section, we discuss the behavior of the response in the time Sector a < 1, which is governed by rare 
events. 



IX. RESPONSE IN THE TIME SECTOR a(t,t w ,f) < 1 FROM RARE EVENTS 

As explained before, in the time sector a(t,t w ,f) < 1 i[112|l . the effective dynamics governed by the decimation 
procedure gives no contribution, and the response will thus be governed by rare events. A similar situation was 
already found in the RSRG studies on the Sinai model [2|| and and on the out-of-equilibrium dynamics of the random 
field Ising model |28j . 



A. Description of the rare events 



For the trap model considered here, the 'rare events' that are responsible for the response in the time sector 
a(t,t w ,f) < 1, can be describded as follows : the particle, which is assumed to be trapped in a renormalized trap 
r > R(t w ,0) at time t w in the effective dynamics, has actually a small probability to be at t w in a "small" trap, i.e. 
an already decimated trap r < R{t Wl 0), for two reasons : 

(i) when /x is small but finite, the particle can be found at time t w in a trap r < R w with a probability of order /i, 
as explained in [|J , where the corrections in /i with respect to the effective dynamics of the unbiased trap model were 
studied in details. 

(ii) when t w is large but not infinite, there is a small probability that the particle is doing an excursion (|34I38|) at 
time t w . 

In both cases, the particle that happens to be in a "small trap" at t w will respond to the external field in the time 
sector a < 1. 



B. Correction of order fi to the effective dynamics 

In the previous study on the unbiased trap model 0, we have studied in details the first corrections at order /i to 
the effective dynamics. In particular, we have shown that in the asymptotic time regime, the particle can be found 
with a probability of order fi in the biggest trap S contained in the interval ]M_, M + [ between the two renormalized 

traps around the origin at scale R(t w ,f = 0). In particular, we have obtained that the probability ip^ (r) to be at 
t w in a trap of trapping time r was given in the domain t < R w at first order in \x by 



^(r)9(r < R w ) ~ U^j K 2 (J*pJ (164) 

where we have dropped numerical factors of order 1 (we refer the reader to 9] for precise results). In particular, this 
probability presents an essential singularity at the origin 

^M^M^e-V 3 ? (165) 

which means that it is very unlikely for the particle to be trapped at t w in the biggest trap S contained in the interval 
]M—, M+[ if the trapping time t of S is much smaller than R w . 

As a consequence, the integrated probability up to scale R (with R < R w ) 



R r +oo 
R 



dTif>% (t) ~ fi I dzz 2 K 1 (z)K 2 (z) (166) 
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that is of order /i when a ~ 1, will also present an essential singularity for small a = R/R n 



I 



R 1 



dT ^J T ) ~ M"^" 73 (167) 



C. Probability to be doing tin excursion cit £^ 

_Z. Probability to be doing an unsuccessful excursion at t w 

From the discussion on unsuccessfull excursions, the probability for a particle to be trapped in the vicinity of a 
renormalized trap tq at scale i?, with neighbors at distances l± reads when the external held vanishes l|34fl 

P^iro^X^^t^L^ll+L ( 1 68 ) 

tin or o 

After the average over tq with the measure (|14|) and over the lengths i|13|) , we obtain that the probability to be doing 
an unsuccessfull excursion at time t w reads H110(l 



P out (t w ) = P^iro, i+, A_) ~ fiR^ 1 = ^ti+" (169) 

As expected, the probability of these rare events is very small when t w is large since fi < 1. In addition, there is a 
prefactor /x that makes this probability even smaller in the limit ji — ► 0. 

2. Probability to be doing a successfull excursion at t m 

The particle may be in a successfull excursion at t w , if it belongs to a renormalized trap that gets decimated around 
the scale R w . From the discussion on successfull excursions, the diffusion time for a length I ~ R!^ reads when the 
external field vanishes reads (|A22IA24I) 

t diff (f = 0,l)~ l l~R% (170) 

This time window of width (At) w ~ ~ around t w corresponds in RG scale to a window around R w ~ 

of width Ai? w ~ (At)t„ ~ ti + " As a consequence, the probability to be doing a successfull excursion at t w can 

be obtained from the probability to be in a trap r ~ R w (|14|l times the window width just estimated 

P dzff (t w ) ~ q Rw {R W )AR W ~ /it#* (171) 
The probability is again very small as expected, and happens to have exactly the same scale as (|169fl . 

3. Probability to be in a small trap r during an excursion 

For a particle doing an excursion, we are now interested into the probability ipR XC ( T ) to be in a trap of trapping 
time r. We have obtained above that the total probability to be in a small trap r < R w behaves as 

RW d T r R x :(r) ~ fitP ~ llR^- 1 (172) 

Assuming that the dependence of ip% x S, (t) in r is given by rq(r), i.e. the initial distribution © weighed by the trapping 
time r (i.e. the particle spends in each trap a time proportional to its trapping time), we obtain the estimation 

r£(r)-Rl^- 1) -^e(r<R w ) (173) 
■j-h- 

which is rather flat for small /i — > 0, in contrast with l|165|) . In particular, the integrated probability up to scale R 
(with R < R w ) reads 

f R 

/ dT^(T) ~ R 2 w ^- 1) R 1 -> i = ^^a 1- " (174) 
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D. Contribution to the response of these rare events 



When the initial condition at t w is a 'small' trap of trapping time r, the effective dynamics will become active again 
when R(f, t — t w ) reaches r, i.e. when a reaches < 1, and the corresponding contribution to the response reads 



(<T>) T {t,t w J) = 9{R > t)RX [F = pfRP) 



(175) 



where the condition 9(R > r) means that the trap has been decimated at scale R, and thus the response is given 
in terms of the scaling function (|81|l found before for the case t w = 0. Averaging over r with the total probability 

i ) Rm{' r ) — , 4'r\ t ) + 4 , 'r c ( t ) coming from the two kinds of rare events describded above, we obtain the leading term 
of the response in the sector a — R/R w < 1 as 



(< x >) r 



»(*,*»,/)= / dTi> Rw ( T )(<Y>) r (t,t w J) = RXo(F = pfRt t ) J 



(176) 



Using the estimations of the integrated probabilities found before (|167I174(I . we finally obtain that the contribution 
of excursions dominate for small a < a w , whereas the contribution of the corrections in fi to the effective dynamics 
are dominant for a w < a < 1. The scale of the crossover value a w can be estimated from the equality between l|167|) 
and \ 174(1 at leading order 



1 



(In R u 



(177) 



In conclusion, before the response of the effective dynamics in the sector a > 1 given in (|158I159I162(I . there exists 
a response in the sector a < 1 as a consequence of rare events. For very small a < a w , the response comes from 
the particles doing excursions in anomalously small traps at t w , and it is reduced by a very small prefactor of order 
/xt^ M_1 ^^ +1 ^(t — £ U) )( 1 ~' i )/( /1+1 ) (|174l) . On the other hand, for a w < a < 1, the response is governed by particles 
which are "in delay" with respect to the effective dynamics, as a consequence of \i being finite, and the response is 
reduced by only a factor of /x. 



X. DISCUSSION OF THE FLUCTUATION-DISSIPATION RELATION 
Linear response regime in a given sample in the effective dynamics time sector 



From the diffusion front l|122fl in a given sample characterized by (A + , A_, A + , A_), we obtain at lowest order in F 
the following results for the rescaled mean position 



< x > (t,t w ,f) 



£{t-t w ,f) 
and the rescaled mean sqare displacement 
< (x - x w ) 2 > (t,t w ,f) 

e(t-t w j) 



<Y>(a,F) = 



A+A- 



X + X w + A+A 
2a 



F + 0(F 2 ) 



=< Y 2 > (a,F) = 



\+\- + 



A+A- + A+A 



a 



0(F) 



(178) 



(179) 



So the Fluctuation-Dissipation relation or Einstein relation is valid in the whole time sector a > 1 as long as the 
linear response is valid and reads in unrescaled quantities with F = (3f£{t — t w ,f), 



0f 

< x > (t, t w ,f) ~ ^- < (x - x w ) 2 > (t, t w ,Q) 



(180) 



This is in agreement with the scaling arguments and numerical simulations presented in |22| and with the non-linear 
Fluctuation Theorem discussed in [2Jj |, that proves that the FDT relation is valid in any given sample for arbitrary //. 

The validity of the FDT relation for the trap model in its aging phase is nevertheless quite remarkable, since the 
dynamics is completely out-of-equilibrium : indeed, in the effective dynamics, the weights of the two important traps 
are not given by Boltzmann factors, they don't even depend on the energies of these two traps, but they are given by 
the probabilities to reach one before the other one, and they thus only depend on the distances to the origin! This 
example, with the explicit expressions in a given sample (|178I179(I . thus shows that the validity of the FDT relation 
in the linear response regime does not imply that the system is at equilibrium or even near equilibrium. 
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B. Non-linear response in the asymptotic aging sector a — > oo 



In the asymptotic aging sector a — * oo, which also corresponds to the case where the external field is applied from 
the very begining t w — 0, we have found a very simple relation (|88|l between the scaling functions for the mean 
displacement and the thermal width which is valid for arbitrary F and in particular in the whole non-linear response 
regime. However, this relation found for disorder averaged quantities does not seem to have a simple interpretation, 
since in a given sample, there is not such a relation between the rescaled mean position and the rescaled thermal 
width that are given by (|74|l 

^(t J) l_ e -F(A + + A-) ^ 

<Ax 2 (t,/)> e -^-(l-e-^ + )(l-e-^-) 

e(tj) + ^ ) — (1 _ e -F(A ++ A- ))2 — ( 182 ) 

Moreover, we have not found an equivalent relation when a is finite (jl53|) . Nevertheless, after the average over the 
samples, we have obtained the following simple property at very long times in the aging regime a — ► oo (|114|l . i.e. in 
the non-linear response regime 



lim <*>(*.*"•/) =/3/ (183) 
^°°\<Ax 2 (t,t w ,f)>) 

whereas, for comparison, there is a factor (1/2) in the FDT relation of the linear response regime l|18U|) 



< x > (t,t w J) \ _ (3f 



< Ax 2 (t,t w ,f) > I /-o 2 



As a comparison, in the pure trap model, it is immediate to obtain the equations for the mean displacement 
(d < n >p U re /dt — 2 sinh 4r and for the thermal width d < An 2 > pur e /dt — 2 cosh 4r so that their ratio is simply 



<n> pure (t) _fif 



< An 2 > pure (f) 2 2 

for arbitrary time, in the regime of small asymetry we are interested in (|1(JH 



tanh — ~— (185) 



XI. CONCLUSION 



We have studied in details the dynamics of the one dimensional disordered trap model when an external force is 
applied from the very beginning at t = 0, or only after a waiting time t w , in the linear as well as in the non-linear 
response regime. Using a disorder-dependent real-space renormalization procedure that becomes exact in the limit 
of strong disorder fi — > 0, we have shown that the diffusion front in each sample consists in two delta peaks, which 
are completely out of equilibrium with each other, since their weights represent the probabilities to reach one before 
the other one. The statistics of the positions and weights of these delta peaks over the samples was then used to 
obtain explicit results for many observables, such as the diffusion front, the mean position, the thermal width, the 
localization parameters and the two-particle correlation function. Since the renormalization procedure is defined 
sample by sample, our approach provides a very clear insight into the important dynamical processes. 

From a more general perspective, it seems that up to now, the studies on the response of aging systems to an external 
field have been mainly restricted to the linear response regime |l7l Il8| , which holds for fixed times (t w , t) in the limit 
of vanishing field / — > 0. However, as in the trap model discussed here, it should be expected in a broad class of aging 
systems that, for a fixed small field /, the validity of the linear response regime is limited in the time sector for (t w ,t) 
by a characteristic time t(/) depending on the external field. Indeed, it seems rather natural that an external field, 
even if it is arbitrarily small, will, for sufficiently long times, drive the system into a configurational landscape which 
is completely different from the initial one. So in the asymptotic time regime beyond the characteristic scale r(/), 
the response will always be governed by non-linear effects. For the special case t w — 0, the full response including 
these non-linear effects has already been studied for the Sinai model |2(j| as well as in the coarsening dynamics of 
the random field Ising model [2^|, via the RSRG approach : in both cases, as in the trap model, the field introduces 
a characteristic time separating the linear response regime from a non-trivial aging regime with non-linear effects. 
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This scenario should more generally apply to other coarsening dynamics. However, in numerical studies on domain 
growth processes, to get better results for the linear response regime at large times, it is usually the response to a 
random field that is measured an d n °t the response to a constant field, that would favor one of the phase 

and induces non-linear effects rapidly. Since this choice of a random field for domain growth processes, is in some 
sence the equivalent of a constant field for spin-glasses [3(| , it seems that the non- linear response of spin-glasses will 
be very different from coarsening systems. For instance, in the dynamics of the spherical Shcrrington-Kirkpatrick 
spin-glass model |38j . the magnetic field introduces a characteristic time that separates the aging dynamics of the 
linear response regime from an equilibrium dynamics at large times : here, the magnetic field does not lead at large 
times to a non-trivial aging regime with non-linear effects in the field, but rather gives rise to an interrupted a ging 
phenomenon. However, this scenario is not expected to hold for other spin-glasses such as the usual SK model |38| . 
if one considers the number of metastables states in a field |39| . In conclusion, the understanding of the non- linear 
effects thar arise at large times in the response of aging systems is still very incomplete and should give rise to further 
studies in the future. 
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APPENDIX A: STATISTICAL PROPERTIES OF EXCURSIONS IN THE PRESENCE OF A FIELD 

As explained in the text, to study the excursions in the renormalized landscape in the presence of a field, we have 
to study the following standard problem : what is the probability distribution Pt (x) of the time t of the first-passage 
at x = without having touched the other boundary x = I before, for a pure random walk starting at x with the 
asymetry h = h(f) > © ? 

For i = l, I — 1, the probability distribution Pt(x) satisfies the equation 

dtPt(x) = (1 + h)P x -i(t) + (1 - h)P x+1 {t) - 2P t {x) (Al) 
with the boundary conditions Po(t) = 5(t) and Pi(t) — 0. So the Laplace transform with respect to t 

p + CC 

P x {s) = / dte' st P t (x) (A2) 
Jo 

satisfies 

(1 + h)P x+1 (s) + (1 - /i)P x -i(s) - (2 + s)P x (s) = (A3) 

for x = 1, I — 1 with the boundary conditions Pq(s) = 1 and Pi(s) = 0. 
The solution reads 

Ms) = M^zMm ( A4 ) 
p l + (s) - p l _{s) 

in terms of the roots 



, , 2 + s ± Vs 2 + 4s + 4/j 2 

^ ±(s) = WTh) (A5) 



The series expansion in s then yields the first moments 

,k 



k (x) = / dtt k P t (x) (A6) 







1. Escape probabilities 

For s — 0, the roots become (we assume h > 0) in terms of the bias / (JSJ) 



p+(0) = 1 (A7) 
p_(0) =i^|= e -W (A8) 
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and the probability to reach before I when starting at x thus reads 

-0fx _ -pfl 

6 (x) = P x (s = 0)= - _ m ( A9) 



a. Escape probability along the drift 

When starting at x = 1, the probability to reach x — I without any visit to x — reads 

1 - e-Pf 

p+(lJ) = l-9 Q (x = l) = T — m (A10) 
Here we are interested in the regime (3f <C 1 and I ^> 1, where the escape probability takes the scaling form 

p+(l,f)~jE+(u = Pfl) (All) 



with the scaling function 



E + (u) = (A12) 

1 — e u 



In particular, E + (u — > 0) — > 1 corresponds to the unbiased case where the escape probability is simply l/l. In the 
other limit where u =— * +oo, we have E + (u) ~ u and the escape probability becomes (/3/) i.e. it is proportional to 
the drift and independent of I. 

b. Escape probability against the drift 
When starting at x = I — 1, the probability to escape to without without any visit to x = I reads 

p-(l,f)=9 (x = l-l)= { - 1 - > e _ pfl (A13) 

which varries as it should between pj(l,f — * 0) = l/l for the unbiased case and p^(l,f — > 1) = for the directed 
case. In the regime (3f -€i 1 and I ^> 1, the escape probability takes the scaling form 

p-{l,f)~- l E-{u = [ifr) (A14) 



with the scaling function 



In particular, E (u — > 0) — > 1 corresponds to the unbiased case where the escape probability is simply l/l. In 
the other limit where u — > +oo, we have E~{u) ~ ue~ u and the escape probability becomes exponentially small 
p-(l,f)~/3fe-Pf l . 

1. Mean time for unsuccessful! excursions 



The expansion at first order in s of 1A4(I yields i|A6() the mean time to reach without any visit to x = I when 
starting at x 



26 

a. Unsuccessful Excursion along the drift 

For x = 1, in the limit /3/<l and I 3> 1, the mean time of unsuccessful excursions along the drift takes the scaling 
form 

*£,(/> = 6 1 (x = l)~ lQ(u = 0fl) (A17) 

where the scaling function 

^/ x 1 -2ue- u -e- 2u , A1 . 

interpolates between 

9( M ) = i-^ + 0( M 4 ) (Al9) 
for the unbiased case, where the mean time is (Z/3), and 



6(«) ~ - (A20) 

u— >oo u 

where the mean time is l /(/?/). 

b. Unsuccessful Excursion against the drift 

The unsuccessful excursions against the drift have the same properties 

0us(f,l) = 0t s (f,l) (A2l) 

3. Mean time for the successful excursions 

Similarly, we find that the mean time needed to reach x = I when starting at x = l for a random walk conditioned 
not to visit x — takes the scaling form 

tj iff (f,l)^l 2 D(u = 0fl)) (A22) 

where the scaling function 

u - 2 + (u + 2)e~ u 
u 2 {\ - e~ u ) 

interpolates between 

D(u) = l + 0(u 2 ) (A24) 



D(u) = U ~.Z-^: (A23) 



for the unbiased case, where the mean diffusion time is (I 2 / '6), and 

D(u) = l + 0^) (A25) 

where the mean diffusion time is l/(/3f). 

Similarly, we find that the mean time for a successful excursion against the drift has the same properties 

**//(/> =*£//(/> (A26) 
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